Matrix Coordinate Bethe Ansatz: 
Applications to XXZ and ASEP models 

N. Crampe^fl, E. Ragoucyfl and D. Simon'^@ 

Universite Montpellier 2, Laboratoire Charles Coulomb UMR 5221, F-34095 

Montpellier, France 

^ CNRS, Laboratoire Charles Coulomb, UMR 5221, F-34095 Montpellier 

Laboratoire de Physique Theorique LAPTH 

CNRS and Universite de Savoie. 
9 chemin de Bellevue, BP 110, F- 74941 Annecy-le-Vieux Cedex, France. 

LPMA, Universite Pierre et Marie Curie, 
Case Courrier 188, 4 place Jussieu, 75252 Paris Cedex 05, France 



Abstract 

We present the construction of the full set of eigenvectors of the open ASEP and 
XXZ models with special constraints on the boundaries. The method combines both 
recent constructions of coordinate Bethe Ansatz and the old method of matrix Ansatz 
specific to the ASEP. This "matrix coordinate Bethe Ansatz" can be viewed as a non- 
commutative coordinate Bethe Ansatz, the non-commutative part being related to the 
algebra appearing in the matrix Ansatz. 
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1 Introduction 



Questions of integrability of models in statistical and quantum mechanics have been much 
more studied for periodic systems than for open systems, for which the numbers of particles 
and excitations may vary. However, open boundary conditions have become central in non- 
equilibrium physics, for which exactly solvable models are needed to explore new features. Fast 
advances have been made in the eighties to generalize Yang-Baxter equations to the open case 
(see [1]). However, although the boundaries which preserve the integrability have been classified 
quite easily [21[3] , the computation of the eigenvalues and the eigenvectors for the non-diagonal 
boundaries is a tricky problem. Indeed, no standard method exists to diagonalize completely 
Sklyanin-type transfer matrices with non-diagonal boundary matrices. Recently, progresses 
have been made in this direction |^]-[T0]. but the problem is far from been solved. The present 
paper addresses this question for a simple model that is used both in statistical mechanics and 
in spin chains. 

The Asymmetric Simple Exclusion Process (ASEP) is the simplest model of transport of 
hard-core particles along a one-dimensional lattice, that exhibits non-trivial behaviours. Each 
site of the lattice may be occupied by one particle or empty. On the boundaries, two reservoirs 
add or remove particles with different rates and create a non-zero flux particle from one to the 
other. For a lattice of L sites, the state space is 2^-dimensional and it can be mapped exactly 
to a system of L spin 1/2, the role of the Markov transition matrix being held by the XXZ spin 
chain Hamiltonian. The role of the reservoirs corresponds in this case to non-diagonal magnetic 
fields on the boundaries. Although all the coefficients have very different interpretations and 
may be complex or real depending on the model, the mathematical structure remains the 
same and integrability techniques can be applied in both cases. In the present paper, ASEP 
notations are chosen, without any loss of generality if one assumes that the coefficients may 
become complex. The detailed model is described in section [21 as well as previous integrability 
results on this model. 

Although the models are identical, there exists at least one method that has been devel- 
oped and applied only to the ASEP case: the Matrix Ansatz (or DEHP method) [11]. The 
reason for this fact is that it requires the a priori knowledge of one eigenvalue to build the 
corresponding eigenvector; it is precisely the case for stochastic models for which the existence 
of invariant measures is known. Attempts to establish parallels between Bethe and Matrix 
Ansatze produced various partial results [12l|13] but never a full understanding of the exact 
relation between both approaches. The present paper shows how to build coordinate Bethe 
Ansatz eigenvectors from a Matrix Ansatz-type vacuum state in section [31 

The plan of the paper is the following. We start in section |2] with a brief summary on 
open ASEP models. Then, in section [3l we present our method, that can be viewed as a 
non-commutative coordinate Bethe Ansatz. The non-commutative part is based on the algebra 
appearing in the matrix Ansatz. We conclude in section HI Appendices |A] and [B] gather 
technical results needed in our construction. Finally, appendix O is devoted to the study of 
finite dimensional representations of the matrix Ansatz algebra. 
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2 ASEP Hamiltonians and constraints on boundaries 



The Markov transition matrix for the open ASEP model is given by 

L-l 

i=i 

where the indices indicate the spaces on which the following matrices act non trivially 

It is well-established [TH4T6] that this ASEP model is related by a similarity transformation 
to the so-called integrable open XXZ model [1] . Some care must be taken with the properties of 
the matrix W . In the case of the ASEP with s = 0, the matrix is stochastic but not hermitian. 
In the context of the ASEP, s is the parameter of the generating function of the current of 
particles and, for s 7^ 0, the matrix W is even not stochastic. In the general case encompassing 
both ASEP and XXZ, there is no stochastic nor self-adjoint property for W and left and right 
eigenvectors are generically different. 

Among the whole set of integrable boundaries, there is one very special subset of constraints 
for which the resolution is simplified. This subset, in the notation of the ASEP model, is 
determined by the following constraints between the boundary parameters a, /3, 7, 5 and s but 
also of the bulk parameters p and q [9|[T0|[T6] 

for e, e' = ± and n G {0, 1, 2, . . . , L - 1} (2.3) 
Cj^{x,y)=y/x and c_(x,?/) = — 1. (2.4) 

The four possible choices for the signs e, e' and the corresponding constraint are summarized 
in table [H as well as the corresponding notation used for the Markov matrix W when the 
constraints are satisfied. 
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Table 1: Possible values for the parameters and the constraints imposed by (12. 3p where n takes 
integer values between and L — 1. 

Similarly to W , we denote by K^'^' and K^"^' the boundary matrices K and K when the 
corresponding constraint (12. Sp is satisfied. In the ASEP, the parameters a, /3, 7 and 5 are 
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positive. Then, only the first two constraints in the table [T] may be considered as relevant. 
However, we treat also the last two sets since they become relevant when we map the ASEP 
problem to the XXZ model. 

By numerical investigations in [T7] , it has been established that the whole spectrum of W^'^ 
is given by two different types of Bethe equations. In previous papers [Hllin] , we succeeded in 
computing the eigenvalues and the left or right eigenvectors corresponding to only one type 
of Bethe equation via a generalization of the coordinate Bethe Ansatz [18] . It seems that this 
generalized coordinate Bethe Ansatz is not enough to obtain the missing cases. An advantage 
of this method consists in giving an interpretation of the number n entering in the constraint: 
it is the maximal number of pseudo-excitations in the Ansatz. 

For W£Zi (i-e. = 1), a second method, now called matrix Ansatz (or DEHP method), 
has been developed in [11] to find the second part of the spectrum. We give the outlines of the 
historical method in subsection 13.11 Then, we present the new results in subsection 13.31 which 
consists in a generalization of the matrix Ansatz. It allows us to obtain the second part of the 
spectrum for any W^"" , i.e. for any subset of constraints. 

In [Smnj, we proved also that for a new type of constraints the generalized coordinate Bethe 
Ansatz also provides one part of the spectrum 

c*(a,7)c*,(/3,5) = (^^^ for e, e' = ± and n e {0, 1, . . . , L - 1} (2.5) 
p — q + v — u±. \/{p — q + V — uY + Auv 



2u 



2.6) 



The four possible choices for the signs e, e' and the corresponding constraints are summarized 
in table |2l We denote by W^'^ the matrix W when one constraint (12.51) is satisfied. 



e 


e' 


Constraints 


w 


+ 


+ 




(a,7)c;(/3,5) = e^ 




W++ 








(a,7)c*_(/3,5) = e^ 


- 1 

J 




+ 




r* 


(a,7)c*_(/3,5) = e^ 


1 - 1 

v<? / 


>v+- 




+ 


c*_ 


(a,7)c;(/3,5) = e^ 


I p ] 


>v~+ 



Table 2: Possible values for the parameters and the constraints imposed by (12. 5p . where n is 
an integer between and L — 1. 



As for the matrix W^'^' , in [TU], we computed one part of the spectrum for W^'*^' Markov 
matrices thanks to the generalized coordinate Bethe Ansatz. We believe that a very similar pro- 
cedure as the one presented in following section 13.31 should hold in this case when we transform 
the matrix as in [10]. We will not give in detail the computations in the present paper. 
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3 Matrix coordinate Bethe Ansatz 



3.1 Matrix Ansatz without excitation 

In this subsection, we present the method introduced in [11] to get the eigenvector with van- 
ishing eigenvalue (i.e. the steady state or the invariant measure) for the Markov process with 
s = 0. This state corresponds exactly to the part of eigenspace we do not obtain from the 
coordinate Bethe Ansatz for W£Zi in [9|[T0]. 

A particular care is required since it involves various operators acting on different vector 
spaces. The physical state space, written all along the paper Ti, is 2^-dimensional and the 
canonical basis can be indexed by the occupation numbers G {0, 1} (resp. spin Si E { — |, |} 
in the XXZ language) on each site i. All the vectors of "H are written in the ket notation |-) 
and all the dual vectors are written in the bra form (-1. 

The matrix Ansatz states that the steady state 1$) of the ASEP with s = has components 
given by: 

{nT2 ...tl\<I>) = {{V,\ H {t,D + (1 - T,)E) \V2)) , (3.1) 

l<i<L 

where the arrow means that the product has to be build from left to right when the index 
i increases. One has for example (011001|$) = {{Vi\EDDEED\V2)) for L = 6. The non- 
commuting matrices D and E act on an abstract auxiliary vector space V, that is different 
from "H. The vector |Vi)) lies in this space V (as any vector written with a ket notation |-))), 
whereas the vector is in its dual V* (as any vector written with a bra notation ((-1). 

Eq. fl3.ip expresses the components of 1$) G "H as a scalar product between vectors on the 
abstract auxiliary space V. 

One checks that |$) is an eigenvector of W£Zi with vanishing eigenvalue if the two matrices 
D et E acting on V and the two boundary vectors satisfy the commutation rules: 

qED - pDE = D + E, (3.2a) 
{{Vi\{aE - -fD + 1) = , (3.2b) 
{/3D - 6E + 1)\V2)) = . (3.2c) 

These three relations reduce quadratic relations in D and E to linear expressions in D and E 
in the bulk and linear relations in D and E to scalar expressions on the boundaries; it allows 
one to determine recursively all the components of the eigenvector |$) and it does not need 
an explicit representation of the algebra. An irrelevant minus sign relatively to the standard 
notation for the matrix Ansatz has been introduced in f l3.2p for later convenience. In particular, 
the matrix Ansatz gives an easy access to correlation functions with standard transfer matrix 
techniques [TT] . 

To prove easily that 1$) is an eigenvalue of Wj^Zi, we firstly rewrite |$) as follows 

1$) = ((^il ( ^ ) ® ( ^ ) ® ■ ■ ■ ® ( ^ ) 1^2)) (3.3) 

The convention used in this notation will be used all through the paper and corresponds to 
the following interpretation: a vector ( ^ ) has operator entries instead of scalar ones (two 
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(3.4) 



dimensional module over the endomorphisms of V), the tensor product of such elements is 
an element of a 2^-dimensional module over the endomorphisms of V, whose components are 
products of size L of operators E and D with the usual tensor rule: 

/aa'\ 
^ a\ ( a!\ ah' 
h)®\h')^ ha' 
\hh'J 

Let us emphasize on the order of the non-commuting operators in the vector of the r.h.s. of 
the previous equation. The right action of IV2)) G V produces a 2^-dimensional vector whose 
components are vectors of V, that is further projected through the left action of ((Vi| G V* to 
a 2^-dimensional vector whose components are complex numbers and is thus identified to "H. 
For example, the following vector 



f{{V,\aa'\V2))\ 
{{Vi\ab'\V2)) 
{{Vi\ha'\V2)) 

\{{Vi\hh'\V,)) / 



(3.5) 



is a 4-component vector with complex entries. 

Secondly, using definitions (12. 2p and equations (13. 2p . we show that 
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-1 



(3.6) 



1^2)) (3.7) 



The first equation encodes four relations on the endomorphisms D and E and each equation on 
the second line two relations between vectors of V or V*. Since this type of result is central in 
this paper and to be pedagogical, we explain in detail the computation for the second relation 
in diZD 



K 



E 
D 



-6 /3 
5 -(3 



E\V2)) 
D\V2)) 



i(3D - SE)\V2)) 
{6E - f3D)\V2)) 



(3.8) 
(3.9) 



where to go from (13. 8 p to (13. 9 p we used (I3.2cp . Let us remark that we denote by 1 the identity 
operator acting on V. 

Finally, remarking that the matrices u, K and K do not act in the auxiliary space V where 
are ((Vi| and IV2)), we obtain telescopic terms that can be simplified to get 



^Lril$) = o 



(3.10) 



The goal of this paper consists in generalizing (13. 3 p to obtain also eigenvalues and eigen- 
vectors for all the VF^'*^ . For that purpose, we add some pseudo-excitations above the previous 
eigenvector as it is done in usual Bethe Ansatz. Before giving the matrix Ansatz with n pseudo- 
excitations, we need to introduce some vectors playing the role of these excitations as well as 
some of their properties. 
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3.2 Pseudo-excitations and properties 

Let us introduce the notation: 

M = ^ (3.11) 

q-p 

and the vectors 

H")>=(Xm) I'W> = (-i) W")> = (p-')(„z) + m) P.12) 

where the notation |.) (ket in bold) stands for 2-component vectors in which the entries are 
operators on V (to differentiate them from vectors |.) G "H introduced in section [XT| which 
contain complex numbers and are obtained after projection with ((Vi| and IV2))). 

If the non-commuting matrices D and E satisfy relation f l3.2ap . we get the following relations 



w\oj{u)) ® \uj{u)) =\uj{u)) (g) \t{u)) - \t{u)) (g) \uj{u)) (3.13a) 

w\u{v)) ® |a;(v-)) =\uj{v)) ® \t{v-)) - \t{v)) ® \uj{v-)) (3.13b) 
P P P 

w\uj{u)) (8 \ijo{v)) = - p\u:{u)) (g) \u:{v)) +p\uj{v)) (g |cj(n-)) 

P 

- \t{u)) (g) \u;{v)) + \u;(u)) (g> \t{v)) (3.13c) 

w\u}{v)) (g \uj{u-)) = - q\uj{v)) ® \u){u-)) + q\u;{u)) ® \oj{v)) 
P P 

+ \u;{v)) ® \t{u-)) - \t{v)) ® \u:{u-)) (3.13d) 
P P 

where u and v are still arbitrary numbers that will be fixed later. The proof of relations f l3.13p is 
straightforward: one projects each relation on the four components then one uses the definition 
of the [u] numbers as well as the commuting relation fl3.2ap . Let us remark that relation ( I3.13ap 
for M = 1 is similar to relation (13. 6p . Relations (I3.13P deal with the bulk part of W. We study 
now these vectors on the boundaries. 

Let us introduce the following functions 

\^{u,v) = and \_^{u,v) = —u — v . (3-14) 

In the following section, we impose the following relation on the right boundary, for e = ±, 

K\uj{u))\V2)) = X,{f3,6)\oj{u))\V2)) - \tiu))\V2)) (3.15) 

where, as explained previously, | V2)) means that IV2)) is right-applied to each entry of 

\oj{u)). This relation provides generally two different constraints {K is a two-by-two matrix) 
on vectors of V except if n = — l/c_e(/3, S). In that case, the remaining constraint becomes 

iPD - 6E + 1)\V2)) = for u = -l/c_,(/3, 5) . (3.16) 

We recover the relation f l3.2cp for both values of Ae(/3, S). 
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Similarly, on the left boundary, we impose the following relation 

K{{V^\\uj{u)) = A,(a,7)((Ki|k(^^)) + m\\t{u)) (3.17) 

where the notation |a;(u)) means that ((Vi| is left-applied to each operator entry of \ijj{u)). 
With a special choice of u, the two constraints reduce again to only one, given by 

{{Vi\(ae'{E-\u])-—{uD + \u]) + l]=Q for m = -e"'c_,(a, 7) . (3.18) 
3.3 Matrix Ansatz with excitations 

Let us define the tensor product (for modules defined on the endomorphisms of V) of vectors 
\ijj{u)) over the sites i to j and write it as 

mu))l = \u:{u))M^))^^l ■ ■ ■ I^H), • (3.19) 

We now fix an integer n and introduce the state with n — m excitations at the ordered positions 
1 < Xm+i < ■ ■ ■ < Xn < L defined as the (C^)'^'^-vector in "H with the projections: 

/ np.(Xm+l-l) + ... + {Xn-l) 

\Xm+l, ■■■,Xn) = [ \ -] (3.20) 



The overall factor \/qJp is introduced only in order to normalize the Bethe roots. The coeffi- 
cients Um and Vm are related through the recursions: 

_ q _ q /o on 

p p 

and the initial coefficients ui and vi are still arbitrary. These vectors correspond to states 
where m excitations have left the system (through the left boundary). To clarify the notation, 
the state with no excitation corresponds to m = n and is given by 

m = {{V,\ 11^2)), (3.22) 

and, for Un+i = 1, we recover the state |$) defined in (13. 3p . 

We need also to introduce some other definitions concerning the set on which we are going 
to sum in our Ansatz. The set Gm is a full set of representatives of the coset BCn/BCm 
(Go = BCn, by convention) and BCm is the Weyl group, generated by transpositions 
(Tj, j = 1, . . . ,m — 1 and the reflection ri (for details, see appendix [X]). It acts on a vector 
k= {ki,...,kn) of C": 

= ( ki, /c2, . . . , kn) , feg-j = (^1, . . . , kj^i, kj, . . . , kn) ■ (3.23) 
We introduce also the following truncated vector kl^\ for < m < n and g G Gm, 

k^r^ = {kgim+l),...,kgrr) . (3.24) 

We are now in position to state the main result of this paper which provides eigenvalues 
and eigenstates of W^^i^n matrix W with constraint c^{a,'y)c^'{/3,6) = . 



Theorem 3.1 The vector 

n 

1^") = E E E 4"^^ ^^'^"^ 1^-+!' • • • '^'^) ' (3-25) 

m=0 Xm+i<---<x„ geGm 

is an eigenstate of W^_i_^ with eigenvalue 

^I'l'i_„ = A_,(a,7) + A_,.(/3,5) + f^A(e^'=0 where A{z) = ^ (^z + - p - q (3.26) 
if the following relations are fulfilled: 

a) The coefficients ui and v\ entering in definition ^3. 2(A) are 

Ml = -e-^c,(a,7) and vi = — [^\ {p-q+5-l3±^J {P - 5 + q - pf + A^S) (3.27) 

h) The non-commuting elements E and D obey 

qED - pDE = E + D (3.28) 
{f3D - 6E + 1)\V2)) = (3.29) 
{{Vi\(ae\E - M) - ^e~\D - [1/u,]) + l) = . (3.30) 



c) The coefficients A^J^^ verify 



= s (e*'=''^ e'^^^f^+i)) , (3.31) 
_ y(m)(gifc,i^ _ _ . ,e^'=^'")A("*-i) , (3.32) 



where 



S{Zi,Z2) = - '^l^^'^^l a(2;^,2;2) = ^ f f ^ J^. \ ^2 - Z1Z2 - 1 ] , (3.33) 

a{z2,zi) Z1Z2 -1 \\ \ p \j q; 



and 



r)(m-i) 2 _ 1 

tM(^^, . . . , ^„) = ^ ,^^_i - , , ^ , , , (3.34) 



Pl{z^)V{{z^) YljJ aiZm, Zj)a{Zj, 1/Zm) 



D^r'^ = (ae%^ + ^-a+p-q-^), (3.35) 



V.^iz) = A{z) + (A^(a,7) + 7)(l - -./^) + (A^(«,7) + «) (l - -./^) , (3.36) 

\ zv qJ \ zv pJ 



piz =2; + ^ . (3.37) 

Jpq vi - U2 
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d) The pseudo-momentum kj must satisfy the following Bethe equations, for I < j < n, 
with 

Vtiz) = A{z) + (A^(/3,5) + /3) (l - - ./^) + (A^(/3,5) + - - ./^) • (3.39) 
Before giving the proof of this theorem in subsection 13. 4[ we make some remarks on the theorem: 



Remark 3.1 The first relation in (13.271) is equivalent, via constraint (12. 3p . to 

n„+i = -l/c,,(/3,5). (3.40) 
The sign in the definition of vi is irrelevant. 

Remark 3.2 The algebra generated by E and D is very closed to the one introduced in [TT]: 
the only difference lies in equation (13.301) . We study the finite dimensional representations of 
this algebra in appendix [Cl which gives intriguing relations with the second set of constraints 

Remark 3.3 A consequence of flX5^ (for m = 1) and A^^^^ = A^g \ is 

^(0) _ T^'Ke'''') ^(0) .3 

This relation with (I3.3ip allow us to express A^g^^ for any g G BCn in terms of Af''^ (where the 
subscript 1 stands for the unit of BCn group). Finally, using recursively (I3.32p . we can express 
all the coefficients A^J^^ in terms of only Af'\ This last coefficient is usually chosen such that 
the eigenfunction be normed. 

Remark 3.4 In our previous work [10], we found via the coordinate Bethe Ansatz the eigen- 
functions of W^'"^ with eigenvalues 

n 

= A,(«,7) + A,,(/3,5) + 5^ A(e^"0 (3-42) 

Remark the change of signs in the index of both A in these eigenvalues (I3.42p in comparison 
with (I3.26p . The pseudo-momentum pj satisfy the following Bethe equations 

n S{e^^\ e^^^)S{e-^^^,e^n = e^^^^^ ,,^1,^^^! wS'^l^ ' J = 1, 2, • • • , n (3.43) 
7=1 Vi (e 'P'Wl (e 'P^) 

By numerical investigations [TT], it has been conjectured that the set of eigenvalues given by 
(I3.42p . (I3.43P with n pseudo-excitations and by (I3.26p . (I3.38P with L — \ — n pseudo-excitations 
give the complete spectrum of W^"" . Therefore, our previous results [10] with the results of this 
paper seem to provide the complete spectrum for W!^'' as well as the associated eigenstates. 
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3.4 Proof of the main theorem 13.11 

To prove the theorem, we show that the following equation holds 

Wrl.^J^n) = Stl-n\^n) ■ (3.44) 

The proof is very similar to the one we performed in [10] for the generalized coordinate Bethe 
Ansatz, except that entries of vectors are now non-commuting operators. It consists 

in projecting equation fl3.44p on the different . . . ,2:„) and to prove that each projection 

is true if the conditions ll])-|d]) of the theorem hold. We write only the projections leading to 
independent relations (one can check that the remaining ones do not lead to new relations). 

On \xxi . . . , Xn) for (xi, . . . , Xn) generic (i.e. 1 < Xi, x„ < L and 1 + Xj < Xj+i) 

Before performing this projection, let us remark that, using relations ( I3.13p . ( I3.15P and 
(I3.17P as well as the conditions |b]) of the theorem, we can show that 

W^2'Ii_„|xi, . . . ,a;„) = (^A_£(a,7) + A_,/(/3,5)j . . . (3.45) 

n 

+ X] ~ 1^1' • • • ) ^n) + \/pg(| . . . , + 1, . . .) + I . . . , Xj— 1, . . .)) 

i=i 

We remind that relations fl3.15p and ( I3.17p are valid only if 



Ml = -e ^Ce(a,7) and m„+i = r^-r- (3.46) 

which is the first relation in[aj) (see also remark WA\ . In addition of that, the recursion relation 
between the m's fl3.2ip implies the constraint Ce(a, 7)ce/(/3, 5) = which is the constraint 

for Wrl^-n- 

Finally, the projection on . . . , Xn) (for generic xi, . . . , x„) of equation fl3.44p holds if the 
energy takes the form (13.260 . 

On |a;i, . . . , Xn) with a3j_|_i = 1 + Xj (and xi, . . . , Xj_i, Xj+2, ■ ■ ■ ,Xn generic) 



Using (I3.13bp and this projection, we get a relation fl3.3ip between A^g^ and A^J^. As ex- 



pected, the expression of the scattering matrix is similar to the periodic case since the bound- 
aries are not involved in this process. 

On \xm+i • • • , Xn) {xm+1, ■ ■ ■,Xn generic and m > 1) 

Before performing this projection, we need to know how the left boundary matrix K acts 
on the vectors \u{vm)) and \u{um)). By direct computation, using relation fl3.30p . we show that 

Kr-l-nm\u^iVrn))=A^r'H{Vl\M^^ (3.47) 

Kr-i-n{mu:{um)) = At''~'\{V,\\u:{Um))+Ct~'\{V^\\u;{vrr._^)) + {{Vi\\t{Vm)) (3.48) 
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where the same notational convention is used as for ( I3.17P and where -D^™" is defined in ( ]3.35p 
and 

Ai™""^^ = (ae'^VmUm+i + (7 + P - Q)vm - aUm+1 - 7e"'*) (3.49) 

^(m~i) ^ (1 - e'uni){avm-i + ^^^^^ 

^m — 1 

^(m-i) ^ (1 - e'u^){au^ + -fe-') 

Using these relations, we get finally that the projection on \xm+i ■ ■ ■ , Xn) provides the following 
constraints, for any g G Gm, 

m 

j^[rn-i) ^(--Dg^fc^M^) + (aM _ A_,(a,7) - 5];A(e^'=-))4™) = 0. (3.52) 

heH„^ j=l 

where Hm = BCm/BCm-i (see Appendix 

We are going to demonstrate that this last constraint fl3.52p holds if relations fl3.3ip and 
(13321) are trutQ. We start by remarking that a consequence of latter equations fl3.3ip and fl3.32p 
is 

{1 l<a<m-l, 
"""gH;.C:a:rr"' g(^"-.^"-") (3.53) 

Then, using again f l3.32p to express now A^f^^ in terms of A^J^"^^ and using f l3.53p to express 
A'-^'^'^ {h e Hm) in terms of A^g'^~^\ relation f l332|) becomes the functional relation 



^m^^)Mz,) ^^^^^ ^^^^^^^^ 1^ ^ z,^^^^^ ^^1 ^^^^^^^^ 

m 

= A_,(a,7)-A("*) + ^A(z,), (3.54) 

i=i 

where Zj stands for exp{ikgj) and the functions are defined in f l3.33p - fl3.36p . To prove this last 
relation f l3.54p . let us introduce the following function 

F(™)(z) = ^^'^'^P'^'^ fla{z,zMzi,-)- (3.55) 

A{z) {2z - - ^) l\ 

Then, one can prove that fl3.54p is equivalent to Xlrcsidues of Ft'") F^"^\z) = (see appendix IB] for 
the complete list of its residues), which finishes to demonstrate that constraint fl3.52p is verified 
if relations f l3.3ip and f l3.32p are true. 



^This demonstration is similar to the one done in }10j but we give here again for completeness of the present 
proof. 
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On |1, Xm+i • • . , Xn) (xm+i, . . .,Xn generic and m > 1) 

This projection provides a second relation between the coefficients from the level m — 1 and 
m. We obtain the following constraint, for any g G Gm, 

m 

(Vwe*"''''" + Ai""'^ - A_,(a,7) - q - ^ A(e*'^«o)4"~'^^''^'"" + ^^^4"^ = -(3.56) 

heH„^ j=l 

We are going to prove that this constraint is satisfied if relations f l3.3ip and fl3.32p following 
similar demonstration as previously. Using the previous projection f l3.52p (already proven) then 
relations f l3.3ip and fl3.32p . we prove that projection (13.560 becomes the following functional 
relation 



^3 



-Vl{zj)'p\{z^) TT a{zj, Zi)a{zi, — ) H ^ ^ TT a( — , Zi)a{zi, zj] 

1 ^ ±± J r- 1 — Zf -^-^ Zi 

1=1 ■> J e=i J 



m m 

= ( J2 A(^.) + 7) - a!"'-'^ + ^) ( E ^(^^O + ^-("' ^) - ^S'^O " C'f™^^!'"-^^ (3.57) 

The function to consider is now (z) = y/pq z F'^'"^\z) . Finally, we prove that functional 
relation (I3.57P is equivalent to Xlrcsiducs of g(™) G^^\z) = (see appendix [B] for the computation 
of the residues). 

On \xi. . . , Xn-i, L) (xi . . . , Xn-i generic) 

To perform this projection, we need to know the action of the right boundary on one 
pseudo-excitation 

Kr\^Ju{vr.))\V,)) = AL|u;K))|y2)) - \t{vr,))\V2)) (3.58) 

where the same convention is used as for (I3.15P and where = 5{vn — 1) and "Wn = ( M Vi 

with Vl given by (I3.27p . 

Then we can prove that this last independent projection holds if the so-called Bethe equa- 
tions (I3.38P are satisfied. 

4 Conclusion 

The previous sections present the construction of the full set of eigenvectors of the ASEP and 
the XXZ spin chain with special constraints on the boundaries. The method combines both 
recent constructions of coordinate Bethe Ansatz for the same set of constraints [9l[T0] and the 
old method of the matrix Ansatz [11] specific to the ASEP. Although computations have been 
showed only one set of special constraints, the construction should be transposed without effort 
to the second set of special constraints discovered in [10] . Left eigenvectors are also very simple 
to build using the same tricks. 
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A first intriguing feature of the matrix coordinate Bethe Ansatz for the first set of special 
constraints fl2.3p is presented in appendix O finite dimensional representations of the matrices 
D and E can be found only if the second set of constraints (12. 5p is satisfied. Up to now, 
however, we do not have simple explanations of this fact. 

The matrix Ansatz has proven to be useful, at least in the case of zero excitations, for the 
study of correlation functions [12] since it is reduced to standard transfer matrix techniques 
of one-dimensional statistical mechanics. The same question in the context of Bethe Ansatz 
is notably difficult and has lead to many different approaches such as the quantum inverse 
scattering method. It would be interesting to investigate whether the present formulation may 
simplify the study of correlation functions, either in the present case with boundaries or in the 
standard periodic case. 

The proofs presented here seem to indicate that the matrix Ansatz state plays the role of 
a new vacuum state, although it is highly non-trivial and do not factorize, in the context of 
the open XXZ spin chain. The standard coordinate Bethe Ansatz approach has been used 
but it would worth knowing if algebraic Bethe Ansatz could be adapted to obtain the same 
eigenvectors, as in 

Numerical evidence tends to show that the set of eigenvectors is now complete and gives a 
synthetic view of the BCn structure of the eigenstates. The constructions for the two sectors 
of the spectrum are similar on their structure but very different from the point of view of 
the reference vacuum state, although a physical interpretation for the ASEP gives some hints 
pUjl^ . One may hope that a further understanding of the passage from one sector to the other 
may allow one to couple both sectors and study ASEP and XXZ spin chains with boundaries 
out of the special constraints that allowed the present framework. 

Finally, the matrix Ansatz approach was found to be useful for various stochastic models 
of particles with different types or impurities [2^123], as well as for the study of tableaux in 
combinatorics |21], even if no Bethe Ansatz approach exists yet for these models. It would be 
interesting to know whether our approach may help for these other models, so that integrable 
system methods can be extended to them. 

A Weyl group BCn and cosets 

The Weyl group BCn is generated by the set {ri, ai, . . . , cr„_i} with the following constraints: 

a'j = l = rl , (Jiricriri = riaiTiCJi , a^aj+iaj = aj+iajaj+i . (A.l) 

The subgroup generated by {cti, . . . , cr^-i} is just the symmetric group. We now consider its 
subgroups generated by {ri, ai, . . . , cXm-i}, m < n, which we identify with BCm- 

For g G BCn, we then define the class gBCm = {gh; h G BC^}, called a left coset. It is 
known that the set of all classes gBCn, which is called BCn/BCm, forms a partition BCn'- "we 
can thus define Cm as a full set of representative of BCn/BCm, such that one has the unique 
decomposition BCn = Ugec™ d^Cm- We set, by convention, Cq ~ BCn and G„ = {I}- 

The action of an element g of Cm on a vector fc*-'"-* = [k^+i, ■ ■ ■ , kn) of M"~™ is given by 
^{rn) _ (^kg(^ni+i), ■ ■ ■ , kg{n))- One checks that this action does not depend on the choice of the 
representative g, such that the action of BCn/BCm is well-defined on M""™' without further 
specifications. This definition is useful because the set {k^^^g G Cm} contains one and only 
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one time the vector {ei^ki-^, . . . , ei^_^ki^_^) for any choice ej = ±, 1 < "ij < n and ij ^ ik- For 
example, {k^g~^^\g G Gn-i} = {{K), {-kn), (kn-i), {-kn-i), • • • , (ki), {-ki)}. 

Finally, we introduce Hm which is a full set of representatives of the coset BCm/ BCm-i 
which may chosen as follows 

riCTi . . . (Tm-20"m-l, (TiriCTi . . . (Tm_20"m-1, • • • , (Jm-l ■ ■ - UiriCJi . . . am~2<^m~l} 



B List of the residues of the function F^^^ and G^'^^ 



We list in this appendix the residues of the function F^"^^ defined by fl3.55p and G^''^\ The 
residues of F^"^^ are 



Res(F(™)(^)) 
Res(F("^)(z)) 



=1 



z=i/2j — 1) Zj Zj 



1 1 ™ 1 



Zj 

1=1 J 



pq - 



g\™ A_,(a,7) + a /p, 

p J p ~ Q y Q 

p\™ A„e(«,7) +7 



Res(F('")(2)) 
Res(F('")(z)) 
Res(F('")(2)) 

Res(F(''")(2)) 

The residues of G^'^''{z) at the point z = zq with zq ^ 0,oo are easy to compute: 



2(g-p) 



g — p y p 



[p- q + 2a- 27)pi(- 



l^fpq 



- E ^(^^) - - ^(2A-e(«, 7) + « + 7) 



Res(G(™n^)) 



Res(F('")(^)) 



Since is not a pole, it remains to compute the residue at infinity: 

m / 

Res(G'("^)(^)) _ = 5Z A(^,)A(2fc) + 2 ( Res(F("^)(^)) 



(B.l) 
(B.2) 

(B.3) 
(B.4) 
(B.5) 

(B.6) 



(B.7) 



pA^^) 2A_,(a, 7) + a + 7 + 
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(B.8) 
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C Representations of the algebra generated by E and D 



In this appendix, we study the finite-dimensional irreducible representations of the algebra 
we used to construct the matrix Ansatz since such representations emerge only for special 
constraints already encountered in the study of the open ASEP or XXZ spin chains. In the 
case without excitation, these types of representation have been studied previously in [T5|l25]. 
We will follow similar proofs for the cases with excitations. 

To be self-contained in this appendix, we remind that the algebra needed in these cases is 
defined by 

qED-pDE = D + E, (C.la) 
i^D - 6E + 1)\V2)) = (C.lb) 

{{V^\(ae'iE - [u]) - ^{uD + M) + l) = , (C.lc) 

where u = — e~*c_e(a, 7) = —{p/q)"/c-e'{f3, 6), e, e' G {±}, c+(x, y) = y/x and c_(x, y) = —1. 

In [25], it is proved that, to satisfy the 'bulk' part ( lC.la|) . the non-commuting elements E 
and D must take the following form (up to similarity transforms) 

where is the dimension of the representation, a is a free parameter and Eij is the elementary 
matrix with 1 in the entry and otherwise. 

Following the arguments of [25], the 'boundary' conditions ( IC.lbp and (IC.lcP imply that 
there exist integers k and i between 1 and A^ such that 

A^ being the size of the irreducible representation, one gets the constraint \k — i\ = N — 1. In 
the case = A^ > 1 and i = 1, one has 

((ri| = (l,0,...,0) and 1^2)) = (0,..., 0,1)* {{V.lD^V^}) = , \f j . (C.6) 

Now, since fIC.laP implies that 

n 

g'^^;^)" = (l+pD)"E + ^g^^i(l+pD)"^^D^' Vn, (C.7) 

e=i 

(1C.6I) shows, together with (IC.lbl) . that all words (in E and D) vanish. Thus, this case must 
be excluded. 
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It remains the case k = 1 and £ = > 1 for which equations (1C.4P and (IC.SP become 

l3a^ + {/3-6 + q-p)a-6 = {C.t 




N-l\ 2 / , / \ 




+ [^^e^-^s+^-p)\-{-] 1-7 = 0. (C.9) 

Then, using definition (I2.6I) . we get, for any choice for u, 

s / \ N-l 

a = cl(f3,6) and — (^) =c!,(a,7) with r, r' = ± . (C.IO) 
a Vgy 

Therefore, a finite dimensional representation exists if there exist two signs r and r' such that 
the following relation is true 

c;(«,7)c;,(A5) = e^(^) (C.ll) 

where the functions have been defined in ( 12. 6p . 

For real parameters and u,v > 0, one has c*_{u, f ) < and c^(m, f ) > (since for uv > one 
has >y {p — q + V — uY + 4ni; > — g + f — n|). Thus, for ASEP models, only the constraints 
corresponding to r = r' have to be considered (as it is the case for the first choice). We recognize 
in the framework of matrix Ansatz the same constraints (I2.5p that appear for coordinate Bethe 
Ansatz, although the relation between both approaches are very different. 
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XXZ models with special constraints on the boundaries. The method combines both 
recent constructions of coordinate Bethe Ansatz and the old method of matrix Ansatz 
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1 Introduction 



Questions of integrability of models in statistical and quantum mechanics have been much 
more studied for periodic systems than for open systems, for which the numbers of particles 
and excitations may vary. However, open boundary conditions have become central in non- 
equilibrium physics, for which exactly solvable models are needed to explore new features. Fast 
advances have been made in the eighties to generalize Yang-Baxter equations to the open case 
(see [1]). However, although the boundaries which preserve the integrability have been classified 
quite easily [21[3] , the computation of the eigenvalues and the eigenvectors for the non-diagonal 
boundaries is a tricky problem. Indeed, no standard method exists to diagonalize completely 
Sklyanin-type transfer matrices with non-diagonal boundary matrices. Recently, progresses 
have been made in this direction |^]-[T0]. but the problem is far from been solved. The present 
paper addresses this question for a simple model that is used both in statistical mechanics and 
in spin chains. 

The Asymmetric Simple Exclusion Process (ASEP) is the simplest model of transport of 
hard-core particles along a one-dimensional lattice, that exhibits non-trivial behaviours. Each 
site of the lattice may be occupied by one particle or empty. On the boundaries, two reservoirs 
add or remove particles with different rates and create a non-zero flux particle from one to the 
other. For a lattice of L sites, the state space is 2^-dimensional and it can be mapped exactly 
to a system of L spin 1/2, the role of the Markov transition matrix being held by the XXZ spin 
chain Hamiltonian. The role of the reservoirs corresponds in this case to non-diagonal magnetic 
fields on the boundaries. Although all the coefficients have very different interpretations and 
may be complex or real depending on the model, the mathematical structure remains the 
same and integrability techniques can be applied in both cases. In the present paper, ASEP 
notations are chosen, without any loss of generality if one assumes that the coefficients may 
become complex. The detailed model is described in section [21 as well as previous integrability 
results on this model. 

Although the models are identical, there exists at least one method that has been developed 
and applied only to the ASEP case: the Matrix Ansatz (or DEHP method) [11]. The reason for 
this fact is that it requires the a priori knowledge of one eigenvalue to build the corresponding 
eigenvector; it is precisely the case for stochastic models for which the existence of invariant mea- 
sures is known. Attempts to establish parallels between Bethe and Matrix Ansatze produced 
various partial results [T2l[T3] (see however [T3] for the case of periodic boundary conditions) 
but never a full understanding of the exact relation between both approaches. The present 
paper shows how to build coordinate Bethe Ansatz eigenvectors from a Matrix Ansatz-type 
vacuum state in section [31 

The plan of the paper is the following. We start in section [2l with a brief summary on 
open ASEP models. Then, in section [3l we present our method, that can be viewed as a 
non-commutative coordinate Bethe Ansatz. The non-commutative part is based on the algebra 
appearing in the matrix Ansatz. We conclude in section [H Appendices \M and [Bl gather 
technical results needed in our construction. Finally, appendix [O is devoted to the study of 
finite dimensional representations of the matrix Ansatz algebra. 



1 



2 ASEP Hamiltonians and constraints on boundaries 



The Markov transition matrix for the open ASEP model is given by 

L-l 

i=i 

where the indices indicate the spaces on which the following matrices act non trivially 

It is well-established [T5HT8] that this ASEP model is related by a similarity transformation 
to the so-called integrable open XXZ model [1] . Some care must be taken with the properties of 
the matrix W . In the case of the ASEP with s = 0, the matrix is stochastic but not hermitian. 
In the context of the ASEP, s is the parameter of the generating function of the current of 
particles and, for s 7^ 0, the matrix W is even not stochastic. In the general case encompassing 
both ASEP and XXZ, there is no stochastic nor self-adjoint property for W and left and right 
eigenvectors are generically different. 

Among the whole set of integrable boundaries, there is one very special subset of constraints 
for which the resolution is simplified. This subset, in the notation of the ASEP model, is 
determined by the following constraints between the boundary parameters a, /3, 7, 5 and s but 
also of the bulk parameters p and q [9|IT0 | [T7 t [T8] 

for e, e' = ± and n G {0, 1, 2, . . . , L - 1} (2.3) 
Cj^{x,y)=y/x and c_(x,?/) = — 1. (2.4) 

The four possible choices for the signs e, e' and the corresponding constraint are summarized 
in table [H as well as the corresponding notation used for the Markov matrix W when the 
constraints are satisfied. 



c,(a,7) 




Constraints 


w 


c+(tt,7) = I 


c+(/3,5) = | 




;?) 


= 1 


' ' n 


C-(«,7) = -1 


c_(/3,5) = -l 


- (?) 


= 1 




c_(a,7) = -1 


c+(/3,5) = | 




(C 


= 1 




c+(a,7) = I 


c_(/3,5) = -l 


7 


(C 


I =1 


n 



Table 1: Possible values for the parameters and the constraints imposed by (12. 3p where n takes 
integer values between and L — 1. 

Similarly to W , we denote by K^'^' and K^"^' the boundary matrices K and K when the 
corresponding constraint (12. Sp is satisfied. In the ASEP, the parameters a, /3, 7 and 5 are 



/ \ 

— g p 

q —p 

loo 00/ 
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positive. Then, only the first two constraints in the table [T] may be considered as relevant. 
However, we treat also the last two sets since they become relevant when we map the ASEP 
problem to the XXZ model. 

By numerical investigations in [in] , it has been established that the whole spectrum of W^'^ 
is given by two different types of Bethe equations. In previous papers [Hllin] , we succeeded in 
computing the eigenvalues and the left or right eigenvectors corresponding to only one type 
of Bethe equation via a generalization of the coordinate Bethe Ansatz |20] . It seems that this 
generalized coordinate Bethe Ansatz is not enough to obtain the missing cases. An advantage 
of this method consists in giving an interpretation of the number n entering in the constraint: 
it is the maximal number of pseudo-excitations in the Ansatz. 

For W£Zi (i-e. = 1), a second method, now called matrix Ansatz (or DEHP method), 
has been developed in [11] to find the second part of the spectrum. We give the outlines of the 
historical method in subsection 13.11 Then, we present the new results in subsection 13.31 which 
consists in a generalization of the matrix Ansatz. It allows us to obtain the second part of the 
spectrum for any W^"" , i.e. for any subset of constraints. 

In [Smnj, we proved also that for a new type of constraints the generalized coordinate Bethe 
Ansatz also provides one part of the spectrum 

c*(a,7)c*,(/3,5) = (^^^ for e, e' = ± and n e {0, 1, . . . , L - 1} (2.5) 
p — q + v — u±. \/{p — q + V — uY + Auv 



2u 



2.6) 



The four possible choices for the signs e, e' and the corresponding constraints are summarized 
in table |2l We denote by W^'^ the matrix W when one constraint (12.51) is satisfied. 
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Constraints 


w 


+ 


+ 




(a,7)c;(/3,5) = e^ 




W++ 








(a,7)c*_(/3,5) = e^ 


- 1 

J 




+ 




r* 


(a,7)c*_(/3,5) = e^ 


1 - 1 

v<? / 


>v+- 




+ 


c*_ 


(a,7)c;(/3,5) = e^ 


I p ] 


>v~+ 



Table 2: Possible values for the parameters and the constraints imposed by (12. 5p . where n is 
an integer between and L — 1. 



As for the matrix W^'^' , in [TU], we computed one part of the spectrum for W^'*^' Markov 
matrices thanks to the generalized coordinate Bethe Ansatz. We believe that a very similar pro- 
cedure as the one presented in following section 13.31 should hold in this case when we transform 
the matrix as in [10]. We will not give in detail the computations in the present paper. 
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3 Matrix coordinate Bethe Ansatz 



3.1 Matrix Ansatz without excitation 

In this subsection, we present the method introduced in [11] to get the eigenvector with van- 
ishing eigenvalue (i.e. the steady state or the invariant measure) for the Markov process with 
s = 0. This state corresponds exactly to the part of eigenspace we do not obtain from the 
coordinate Bethe Ansatz for W£Zi in [9|[T0]. 

A particular care is required since it involves various operators acting on different vector 
spaces. The physical state space, written all along the paper Ti, is 2^-dimensional and the 
canonical basis can be indexed by the occupation numbers G {0, 1} (resp. spin Si E { — |, |} 
in the XXZ language) on each site i. All the vectors of "H are written in the ket notation |-) 
and all the dual vectors are written in the bra form (-1. 

The matrix Ansatz states that the steady state 1$) of the ASEP with s = has components 
given by: 

{nT2 ...tl\<I>) = {{V,\ H {t,D + (1 - T,)E) \V2)) , (3.1) 

l<i<L 

where the arrow means that the product has to be build from left to right when the index 
i increases. One has for example (011001|$) = {{Vi\EDDEED\V2)) for L = 6. The non- 
commuting matrices D and E act on an abstract auxiliary vector space V, that is different 
from "H. The vector |Vi)) lies in this space V (as any vector written with a ket notation |-))), 
whereas the vector is in its dual V* (as any vector written with a bra notation ((-1). 

Eq. fl3.ip expresses the components of 1$) G "H as a scalar product between vectors on the 
abstract auxiliary space V. 

One checks that |$) is an eigenvector of W£Zi with vanishing eigenvalue if the two matrices 
D et E acting on V and the two boundary vectors satisfy the commutation rules: 

qED - pDE = D + E, (3.2a) 
{{Vi\{aE - -fD + 1) = , (3.2b) 
{/3D - 6E + 1)\V2)) = . (3.2c) 

These three relations reduce quadratic relations in D and E to linear expressions in D and E 
in the bulk and linear relations in D and E to scalar expressions on the boundaries; it allows 
one to determine recursively all the components of the eigenvector |$) and it does not need 
an explicit representation of the algebra. An irrelevant minus sign relatively to the standard 
notation for the matrix Ansatz has been introduced in f l3.2p for later convenience. In particular, 
the matrix Ansatz gives an easy access to correlation functions with standard transfer matrix 
techniques [TT] . 

To prove easily that 1$) is an eigenvector of Wj^Zi, we firstly rewrite 1$) as follows 

1$) = ((^il ( ^ ) ® ( ^ ) ® ■ ■ ■ ® ( ^ ) 1^2)) (3.3) 

The convention used in this notation will be used all through the paper and corresponds to 
the following interpretation: a vector ( ^ ) has operator entries instead of scalar ones (two 
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(3.4) 



dimensional module over the endomorphisms of V), the tensor product of such elements is 
an element of a 2^-dimensional module over the endomorphisms of V, whose components are 
products of size L of operators E and D with the usual tensor rule: 

/aa'\ 
^ a\ ( a!\ ah' 
h)®\h')^ ha' 
\hh'J 

Let us emphasize on the order of the non-commuting operators in the vector of the r.h.s. of 
the previous equation. The right action of IV2)) G V produces a 2^-dimensional vector whose 
components are vectors of V, that is further projected through the left action of ((Vi| G V* to 
a 2^-dimensional vector whose components are complex numbers and is thus identified to "H. 
For example, the following vector 



f{{V,\aa'\V2))\ 
{{Vi\ab'\V2)) 
{{Vi\ha'\V2)) 

\{{Vi\hh'\V,)) / 



(3.5) 



is a 4-component vector with complex entries. 

Secondly, using definitions (12. 2p and equations (13. 2p . we show that 



E 
D 



E 
D 



E 
D 



E 
D 



E 
D 



{{Vi 



and K 



E 
D 



-1 



(3.6) 



1^2)) (3.7) 



The first equation encodes four relations on the endomorphisms D and E and each equation on 
the second line two relations between vectors of V or V*. Since this type of result is central in 
this paper and to be pedagogical, we explain in detail the computation for the second relation 
in diZD 



K 



E 
D 



-6 /3 
5 -(3 



E\V2)) 
D\V2)) 



i(3D - SE)\V2)) 
{6E - f3D)\V2)) 



(3.8) 
(3.9) 



where to go from (13. 8 p to (13. 9 p we used (I3.2cp . Let us remark that we denote by 1 the identity 
operator acting on V. 

Finally, remarking that the matrices u, K and K do not act in the auxiliary space V where 
are ((Vi| and IV2)), we obtain telescopic terms that can be simplified to get 



^Lril$) = o 



(3.10) 



The goal of this paper consists in generalizing (13. 3 p to obtain also eigenvalues and eigen- 
vectors for all the VF^'*^ . For that purpose, we add some pseudo-excitations above the previous 
eigenvector as it is done in usual Bethe Ansatz. Before giving the matrix Ansatz with n pseudo- 
excitations, we need to introduce some vectors playing the role of these excitations as well as 
some of their properties. 
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3.2 Pseudo-excitations and properties 

Let us introduce the notation: 

M = ^ (3.11) 

q-p 

and the vectors 

H")>=(Xm) I'W> = (-i) W")> = (p-')(„z) + m) P.12) 

where the notation |.) (ket in bold) stands for 2-component vectors in which the entries are 
operators on V (to differentiate them from vectors |.) G "H introduced in section [XT| which 
contain complex numbers and are obtained after projection with ((Vi| and IV2))). 

If the non-commuting matrices D and E satisfy relation f l3.2ap . we get the following relations 



w\oj{u)) ® \uj{u)) =\uj{u)) (g) \t{u)) - \t{u)) (g) \uj{u)) (3.13a) 

w\u{v)) ® |a;(v-)) =\uj{v)) ® \t{v-)) - \t{v)) ® \uj{v-)) (3.13b) 
P P P 

w\uj{u)) (8 \ijo{v)) = - p\u:{u)) (g) \u:{v)) +p\uj{v)) (g |cj(n-)) 

P 

- \t{u)) (g) \u;{v)) + \u;(u)) (g> \t{v)) (3.13c) 

w\u}{v)) (g \uj{u-)) = - q\uj{v)) ® \u){u-)) + q\u;{u)) ® \oj{v)) 
P P 

+ \u;{v)) ® \t{u-)) - \t{v)) ® \u:{u-)) (3.13d) 
P P 

where u and v are still arbitrary numbers that will be fixed later. The proof of relations f l3.13p is 
straightforward: one projects each relation on the four components then one uses the definition 
of the [u] numbers as well as the commuting relation fl3.2ap . Let us remark that relation ( I3.13ap 
for M = 1 is similar to relation (13. 6p . Relations (I3.13P deal with the bulk part of W. We study 
now these vectors on the boundaries. 

Let us introduce the following functions 

\^{u,v) = and \_^{u,v) = —u — v . (3-14) 

In the following section, we impose the following relation on the right boundary, for e = ±, 

K\uj{u))\V2)) = X,{f3,6)\oj{u))\V2)) - \tiu))\V2)) (3.15) 

where, as explained previously, | V2)) means that IV2)) is right-applied to each entry of 

\oj{u)). This relation provides generally two different constraints {K is a two-by-two matrix) 
on vectors of V except if n = — l/c_e(/3, S). In that case, the remaining constraint becomes 

iPD - 6E + 1)\V2)) = for u = -l/c_,(/3, 5) . (3.16) 

We recover the relation f l3.2cp for both values of Ae(/3, S). 
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Similarly, on the left boundary, we impose the following relation 

K{{V^\\uj{u)) = A,(a,7)((Ki|k(^^)) + m\\t{u)) (3.17) 

where the notation |a;(u)) means that ((Vi| is left-applied to each operator entry of \ijj{u)). 
With a special choice of u, the two constraints reduce again to only one, given by 

{{Vi\(ae'{E-\u])-—{uD + \u]) + l]=Q for m = -e"'c_,(a, 7) . (3.18) 
3.3 Matrix Ansatz with excitations 

Let us define the tensor product (for modules defined on the endomorphisms of V) of vectors 
\ijj{u)) over the sites i to j and write it as 

mu))l = \u:{u))M^))^^l ■ ■ ■ I^H), • (3.19) 

We now fix an integer n and introduce the state with n — m excitations at the ordered positions 
1 < Xm+i < ■ ■ ■ < Xn < L defined as the (C^)'^'^-vector in "H with the projections: 

/ np.(Xm+l-l) + ... + {Xn-l) 

\Xm+l, ■■■,Xn) = [ \ -] (3.20) 



The overall factor \/qJp is introduced only in order to normalize the Bethe roots. The coeffi- 
cients Um and Vm are related through the recursions: 

_ q _ q /o on 

p p 

and the initial coefficients ui and vi are still arbitrary. These vectors correspond to states 
where m excitations have left the system (through the left boundary). To clarify the notation, 
the state with no excitation corresponds to m = n and is given by 

m = {{V,\ 11^2)), (3.22) 

and, for Un+i = 1, we recover the state |$) defined in (13. 3p . 

We need also to introduce some other definitions concerning the set on which we are going 
to sum in our Ansatz. The set Gm is a full set of representatives of the coset BCn/BCm 
(Go = BCn, by convention) and BCm is the Weyl group, generated by transpositions 
(Tj, j = 1, . . . ,m — 1 and the reflection ri (for details, see appendix [X]). It acts on a vector 
k= {ki,...,kn) of C": 

= ( ki, /c2, . . . , kn) , feg-j = (^1, . . . , kj^i, kj, . . . , kn) ■ (3.23) 
We introduce also the following truncated vector kl^\ for < m < n and g G Gm, 

k^r^ = {kgim+l),...,kgrr) . (3.24) 

We are now in position to state the main result of this paper which provides eigenvalues 
and eigenstates of W^^i^n matrix W with constraint c^{a,'y)c^'{/3,6) = . 



Theorem 3.1 The vector 

n 

1^") = E E E 4"^^ ^^'^"^ 1^-+!' • • • '^'^) ' (3-25) 

m=0 Xm+i<---<x„ geGm 

is an eigenstate of W^_i_^ with eigenvalue 

^I'l'i_„ = A_,(a,7) + A_,.(/3,5) + f^A(e^'=0 where A{z) = ^ (^z + - p - q (3.26) 
if the following relations are fulfilled: 

a) The coefficients ui and v\ entering in definition ^3. 2(A) are 

Ml = -e-^c,(a,7) and vi = — [^\ {p-q+5-l3±^J {P - 5 + q - pf + A^S) (3.27) 

h) The non-commuting elements E and D obey 

qED - pDE = E + D (3.28) 
{f3D - 6E + 1)\V2)) = (3.29) 
{{Vi\(ae\E - M) - ^e~\D - [1/u,]) + l) = . (3.30) 



c) The coefficients A^J^^ verify 



= s (e*'=''^ e'^^^f^+i)) , (3.31) 
_ y(m)(gifc,i^ _ _ . ,e^'=^'")A("*-i) , (3.32) 



where 



S{Zi,Z2) = - '^l^^'^^l a(2;^,2;2) = ^ f f ^ J^. \ ^2 - Z1Z2 - 1 ] , (3.33) 

a{z2,zi) Z1Z2 -1 \\ \ p \j q; 



and 



r)(m-i) 2 _ 1 

tM(^^, . . . , ^„) = ^ ,^^_i - , , ^ , , , (3.34) 



Pl{z^)V{{z^) YljJ aiZm, Zj)a{Zj, 1/Zm) 



D^r'^ = (ae%^ + ^-a+p-q-^), (3.35) 



V.^iz) = A{z) + (A^(a,7) + 7)(l - -./^) + (A^(«,7) + «) (l - -./^) , (3.36) 

\ zv qJ \ zv pJ 



piz =2; + ^ . (3.37) 

Jpq vi - U2 
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d) The pseudo-momentum kj must satisfy the following Bethe equations, for I < j < n, 
with 

Vtiz) = A{z) + (A^(/3,5) + /3) (l - - ./^) + (A^(/3,5) + - - ./^) • (3.39) 
Before giving the proof of this theorem in subsection 13. 4[ we make some remarks on the theorem: 



Remark 3.1 The first relation in (13.271) is equivalent, via constraint (12. 3p . to 

n„+i = -l/c,,(/3,5). (3.40) 
The sign in the definition of vi is irrelevant. 

Remark 3.2 The algebra generated by E and D is very closed to the one introduced in [TT]: 
the only difference lies in equation (13.301) . We study the finite dimensional representations of 
this algebra in appendix [Cl which gives intriguing relations with the second set of constraints 

Remark 3.3 A consequence of flX5^ (for m = 1) and A^^^^ = A^g \ is 

^(0) _ T^'Ke'''') ^(0) .3 

This relation with (I3.3ip allow us to express A^g^^ for any g G BCn in terms of Af''^ (where the 
subscript 1 stands for the unit of BCn group). Finally, using recursively (I3.32p . we can express 
all the coefficients A^J^^ in terms of only Af'\ This last coefficient is usually chosen such that 
the eigenfunction be normed. 

Remark 3.4 In our previous work [10], we found via the coordinate Bethe Ansatz the eigen- 
functions of W^'"^ with eigenvalues 

n 

= A,(«,7) + A,,(/3,5) + 5^ A(e^"0 (3-42) 

Remark the change of signs in the index of both A in these eigenvalues (I3.42p in comparison 
with (I3.26p . The pseudo-momentum pj satisfy the following Bethe equations 

n S{e^^\ e^^^)S{e-^^^,e^n = e^^^^^ ,,^1,^^^! wS'^l^ ' J = 1, 2, • • • , n (3.43) 
7=1 Vi (e 'P'Wl (e 'P^) 

By numerical investigations [12], it has been conjectured that the set of eigenvalues given by 
(I3.42p . (I3.43P with n pseudo-excitations and by (I3.26p . (I3.38P with L — \ — n pseudo-excitations 
give the complete spectrum of W^"" . Therefore, our previous results [10] with the results of this 
paper seem to provide the complete spectrum for W!^'' as well as the associated eigenstates. 
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3.4 Proof of the main theorem 13.11 

To prove the theorem, we show that the following equation holds 

Wrl.^J^n) = Stl-n\^n) ■ (3.44) 

The proof is very similar to the one we performed in [10] for the generalized coordinate Bethe 
Ansatz, except that entries of vectors are now non-commuting operators. It consists 

in projecting equation fl3.44p on the different . . . ,2:„) and to prove that each projection 

is true if the conditions ll])-|d]) of the theorem hold. We write only the projections leading to 
independent relations (one can check that the remaining ones do not lead to new relations). 

On \xxi . . . , Xn) for (xi, . . . , Xn) generic (i.e. 1 < Xi, x„ < L and 1 + Xj < Xj+i) 

Before performing this projection, let us remark that, using relations ( I3.13p . ( I3.15P and 
(I3.17P as well as the conditions |b]) of the theorem, we can show that 

W^2'Ii_„|xi, . . . ,a;„) = (^A_£(a,7) + A_,/(/3,5)j . . . (3.45) 

n 

+ X] ~ 1^1' • • • ) ^n) + \/pg(| . . . , + 1, . . .) + I . . . , Xj— 1, . . .)) 

i=i 

We remind that relations fl3.15p and ( I3.17p are valid only if 



Ml = -e ^Ce(a,7) and m„+i = r^-r- (3.46) 

which is the first relation in[aj) (see also remark WA\ . In addition of that, the recursion relation 
between the m's fl3.2ip implies the constraint Ce(a, 7)ce/(/3, 5) = which is the constraint 

for Wrl^-n- 

Finally, the projection on . . . , Xn) (for generic xi, . . . , x„) of equation fl3.44p holds if the 
energy takes the form (13.260 . 

On |a;i, . . . , Xn) with a3j_|_i = 1 + Xj (and xi, . . . , Xj_i, Xj+2, ■ ■ ■ ,Xn generic) 



Using (I3.13bp and this projection, we get a relation fl3.3ip between A^g^ and A^J^. As ex- 



pected, the expression of the scattering matrix is similar to the periodic case since the bound- 
aries are not involved in this process. 

On \xm+i • • • , Xn) {xm+1, ■ ■ ■,Xn generic and m > 1) 

Before performing this projection, we need to know how the left boundary matrix K acts 
on the vectors \u{vm)) and \u{um)). By direct computation, using relation fl3.30p . we show that 

Kr-l-nm\u^iVrn))=A^r'H{Vl\M^^ (3.47) 

Kr-i-n{mu:{um)) = At''~'\{V,\\u:{Um))+Ct~'\{V^\\u;{vrr._^)) + {{Vi\\t{Vm)) (3.48) 
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where the same notational convention is used as for ( I3.17P and where -D^™" is defined in ( ]3.35p 
and 

Ai™""^^ = (ae'^VmUm+i + (7 + P - Q)vm - aUm+1 - 7e"'*) (3.49) 

^(m~i) ^ (1 - e'uni){avm-i + ^^^^^ 

^m — 1 

^(m-i) ^ (1 - e'u^){au^ + -fe-') 

Using these relations, we get finally that the projection on \xm+i ■ ■ ■ , Xn) provides the following 
constraints, for any g G Gm, 

m 

j^[rn-i) ^(--Dg^fc^M^) + (aM _ A_,(a,7) - 5];A(e^'=-))4™) = 0. (3.52) 

heH„^ j=l 

where Hm = BCm/BCm-i (see Appendix 

We are going to demonstrate that this last constraint fl3.52p holds if relations fl3.3ip and 
(13321) are trutQ. We start by remarking that a consequence of latter equations fl3.3ip and fl3.32p 
is 

{1 l<a<m-l, 
"""gH;.C:a:rr"' g(^"-.^"-") (3.53) 

Then, using again f l3.32p to express now A^f^^ in terms of A^J^"^^ and using f l3.53p to express 
A'-^'^'^ {h e Hm) in terms of A^g'^~^\ relation f l332|) becomes the functional relation 



^m^^)Mz,) ^^^^^ ^^^^^^^^ 1^ ^ z,^^^^^ ^^1 ^^^^^^^^ 

m 

= A_,(a,7)-A("*) + ^A(z,), (3.54) 

i=i 

where Zj stands for exp{ikgj) and the functions are defined in f l3.33p - fl3.36p . To prove this last 
relation f l3.54p . let us introduce the following function 

F(™)(z) = ^^'^'^P'^'^ fla{z,zMzi,-)- (3.55) 

A{z) {2z - - ^) l\ 

Then, one can prove that fl3.54p is equivalent to Xlrcsidues of Ft'") F^"^\z) = (see appendix IB] for 
the complete list of its residues), which finishes to demonstrate that constraint fl3.52p is verified 
if relations f l3.3ip and f l3.32p are true. 



^This demonstration is similar to the one done in }10j but we give here again for completeness of the present 
proof. 
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On |1, Xm+i • • . , Xn) (xm+i, . . .,Xn generic and m > 1) 

This projection provides a second relation between the coefficients from the level m — 1 and 
m. We obtain the following constraint, for any g G Gm, 

m 

(Vwe*"''''" + Ai""'^ - A_,(a,7) - q - ^ A(e*'^«o)4"~'^^''^'"" + ^^^4"^ = -(3.56) 

heH„^ j=l 

We are going to prove that this constraint is satisfied if relations f l3.3ip and fl3.32p following 
similar demonstration as previously. Using the previous projection f l3.52p (already proven) then 
relations f l3.3ip and fl3.32p . we prove that projection (13.560 becomes the following functional 
relation 



^3 



-Vl{zj)'p\{z^) TT a{zj, Zi)a{zi, — ) H ^ ^ TT a( — , Zi)a{zi, zj] 

1 ^ ±± J r- 1 — Zf -^-^ Zi 

1=1 ■> J e=i J 



m m 

= ( J2 A(^.) + 7) - a!"'-'^ + ^) ( E ^(^^O + ^-("' ^) - ^S'^O " C'f™^^!'"-^^ (3.57) 

The function to consider is now (z) = y/pq z F'^'"^\z) . Finally, we prove that functional 
relation (I3.57P is equivalent to Xlrcsiducs of g(™) G^^\z) = (see appendix [B] for the computation 
of the residues). 

On \xi. . . , Xn-i, L) (xi . . . , Xn-i generic) 

To perform this projection, we need to know the action of the right boundary on one 
pseudo-excitation 

Kr\^Ju{vr.))\V,)) = AL|u;K))|y2)) - \t{vr,))\V2)) (3.58) 

where the same convention is used as for (I3.15P and where = 5{vn — 1) and "Wn = ( M Vi 

with Vl given by (I3.27p . 

Then we can prove that this last independent projection holds if the so-called Bethe equa- 
tions (I3.38P are satisfied. 

4 Conclusion 

The previous sections present the construction of the full set of eigenvectors of the ASEP and 
the XXZ spin chain with special constraints on the boundaries. The method combines both 
recent constructions of coordinate Bethe Ansatz for the same set of constraints [9l[T0] and the 
old method of the matrix Ansatz [11] specific to the ASEP. Although computations have been 
showed only one set of special constraints, the construction should be transposed without effort 
to the second set of special constraints discovered in [10] . Left eigenvectors are also very simple 
to build using the same tricks. 
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A first intriguing feature of the matrix coordinate Bethe Ansatz for the first set of special 
constraints fl2.3p is presented in appendix O finite dimensional representations of the matrices 
D and E can be found only if the second set of constraints (12. 5p is satisfied. Up to now, 
however, we do not have simple explanations of this fact. 

The matrix Ansatz has proven to be useful, at least in the case of zero excitations, for the 
study of correlation functions [2T] since it is reduced to standard transfer matrix techniques 
of one-dimensional statistical mechanics. The same question in the context of Bethe Ansatz 
is notably difficult and has lead to many different approaches such as the quantum inverse 
scattering method. It would be interesting to investigate whether the present formulation may 
simplify the study of correlation functions, either in the present case with boundaries or in the 
standard periodic case. 

The proofs presented here seem to indicate that the matrix Ansatz state plays the role of 
a new vacuum state, although it is highly non-trivial and do not factorize, in the context of 
the open XXZ spin chain. The standard coordinate Bethe Ansatz approach has been used 
but it would worth knowing if algebraic Bethe Ansatz could be adapted to obtain the same 
eigenvectors, as in 

Numerical evidence tends to show that the set of eigenvectors is now complete and gives a 
synthetic view of the BCn structure of the eigenstates. The constructions for the two sectors 
of the spectrum are similar on their structure but very different from the point of view of 
the reference vacuum state, although a physical interpretation for the ASEP gives some hints 
. One may hope that a further understanding of the passage from one sector to the other 
may allow one to couple both sectors and study ASEP and XXZ spin chains with boundaries 
out of the special constraints that allowed the present framework. 

Finally, the matrix Ansatz approach was found to be useful for various stochastic models 
of particles with different types or impurities [211125], as well as for the study of tableaux in 
combinatorics |2S], even if no Bethe Ansatz approach exists yet for these models. It would be 
interesting to know whether our approach may help for these other models, so that integrable 
system methods can be extended to them. 

A Weyl group BCn and cosets 

The Weyl group BCn is generated by the set {ri, ai, . . . , cr„_i} with the following constraints: 

a'j = l = rl , (Jiricriri = riaiTiCJi , a^aj+iaj = aj+iajaj+i . (A.l) 

The subgroup generated by {cti, . . . , cr^-i} is just the symmetric group. We now consider its 
subgroups generated by {ri, ai, . . . , cXm-i}, m < n, which we identify with BCm- 

For g G BCn, we then define the class gBCm = {gh; h G BC^}, called a left coset. It is 
known that the set of all classes gBCn, which is called BCn/BCm, forms a partition BCn'- "we 
can thus define Cm as a full set of representative of BCn/BCm, such that one has the unique 
decomposition BCn = Ugec™ d^Cm- We set, by convention, Cq ~ BCn and G„ = {I}- 

The action of an element g of Cm on a vector fc*-'"-* = [k^+i, ■ ■ ■ , kn) of M"~™ is given by 
^{rn) _ (^kg(^ni+i), ■ ■ ■ , kg{n))- One checks that this action does not depend on the choice of the 
representative g, such that the action of BCn/BCm is well-defined on M""™' without further 
specifications. This definition is useful because the set {k^^^g G Cm} contains one and only 
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one time the vector {ei^ki-^, . . . , ei^_^ki^_^) for any choice ej = ±, 1 < "ij < n and ij ^ ik- For 
example, {k^g~^^\g G Gn-i} = {{K), {-kn), (kn-i), {-kn-i), • • • , (ki), {-ki)}. 

Finally, we introduce Hm which is a full set of representatives of the coset BCm/ BCm-i 
which may chosen as follows 

riCTi . . . (Tm-20"m-l, (TiriCTi . . . (Tm_20"m-1, • • • , (Jm-l ■ ■ - UiriCJi . . . am~2<^m~l} 



B List of the residues of the function F^^^ and G^'^^ 



We list in this appendix the residues of the function F^"^^ defined by fl3.55p and G^''^\ The 
residues of F^"^^ are 



Res(F(™)(^)) 
Res(F("^)(z)) 



=1 



z=i/2j — 1) Zj Zj 



1 1 ™ 1 



Zj 

1=1 J 



pq - 



g\™ A_,(a,7) + a /p, 

p J p ~ Q y Q 

p\™ A„e(«,7) +7 



Res(F('")(2)) 
Res(F('")(z)) 
Res(F('")(2)) 

Res(F(''")(2)) 

The residues of G^'^''{z) at the point z = zq with zq ^ 0,oo are easy to compute: 



2(g-p) 



g — p y p 



[p- q + 2a- 27)pi(- 



l^fpq 



- E ^(^^) - - ^(2A-e(«, 7) + « + 7) 



Res(G(™n^)) 



Res(F('")(^)) 



Since is not a pole, it remains to compute the residue at infinity: 

m / 

Res(G'("^)(^)) _ = 5Z A(^,)A(2fc) + 2 ( Res(F("^)(^)) 



(B.l) 
(B.2) 

(B.3) 
(B.4) 
(B.5) 

(B.6) 



(B.7) 



pA^^) 2A_,(a, 7) + a + 7 + 
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(B.8) 
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C Representations of the algebra generated by E and D 



In this appendix, we study the finite-dimensional irreducible representations of the algebra 
we used to construct the matrix Ansatz since such representations emerge only for special 
constraints already encountered in the study of the open ASEP or XXZ spin chains. In the 
case without excitation, these types of representation have been studied previously in [T6|l27]. 
We will follow similar proofs for the cases with excitations. 

To be self-contained in this appendix, we remind that the algebra needed in these cases is 
defined by 

qED-pDE = D + E, (C.la) 
i^D - 6E + 1)\V2)) = (C.lb) 

{{V^\(ae'iE - [u]) - ^{uD + M) + l) = , (C.lc) 

where u = — e~*c_e(a, 7) = —{p/q)"/c-e'{f3, 6), e, e' G {±}, c+(x, y) = y/x and c_(x, y) = —1. 

In [27], it is proved that, to satisfy the 'bulk' part ( lC.la|) . the non-commuting elements E 
and D must take the following form (up to similarity transforms) 

where is the dimension of the representation, a is a free parameter and Eij is the elementary 
matrix with 1 in the entry and otherwise. 

Following the arguments of [27], the 'boundary' conditions ( IC.lbp and (IC.lcP imply that 
there exist integers k and i between 1 and A^ such that 

A^ being the size of the irreducible representation, one gets the constraint \k — i\ = N — 1. In 
the case = A^ > 1 and i = 1, one has 

((ri| = (l,0,...,0) and 1^2)) = (0,..., 0,1)* {{V.lD^V^}) = , \f j . (C.6) 

Now, since fIC.laP implies that 

n 

g'^^;^)" = (l+pD)"E + ^g^^i(l+pD)"^^D^' Vn, (C.7) 

e=i 

(1C.6I) shows, together with (IC.lbl) . that all words (in E and D) vanish. Thus, this case must 
be excluded. 
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It remains the case k = 1 and £ = > 1 for which equations (1C.4P and (IC.SP become 

l3a^ + {/3-6 + q-p)a-6 = {C.t 




N-l\ 2 / , / \ 




+ [^^e^-^s+^-p)\-{-] 1-7 = 0. (C.9) 

Then, using definition (I2.6I) . we get, for any choice for u, 

s / \ N-l 

a = cl(f3,6) and — (^) =c!,(a,7) with r, r' = ± . (C.IO) 
a Vgy 

Therefore, a finite dimensional representation exists if there exist two signs r and r' such that 
the following relation is true 

c;(«,7)c;,(A5) = e^(^) (C.ll) 

where the functions have been defined in ( 12. 6p . 

For real parameters and u,v > 0, one has c*_{u, f ) < and c^(m, f ) > (since for uv > one 
has >y {p — q + V — uY + 4ni; > — g + f — n|). Thus, for ASEP models, only the constraints 
corresponding to r = r' have to be considered (as it is the case for the first choice). We recognize 
in the framework of matrix Ansatz the same constraints (I2.5p that appear for coordinate Bethe 
Ansatz, although the relation between both approaches are very different. 
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